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Abstract — We consider error correction procedures designed 
specifically for the amplitude damping channel. We analyze am- 
plitude damping errors in the stabilizer formalism. This analysis 
allows a generalization of the [4, 1] 'approximate' amplitude 
damping code of [1]. We present this generalization as a class of 
[2(M + 1),M] codes for M > 1 and present quantum circuits 
for encoding and recovery operations. We also present a [7, 3] 
amplitude damping code based on the classical Hamming code. 
All of these are stabilizer codes whose encoding and recovery 
operations can be completely described with Clifford group 
operations. Finally, we describe optimization options in which 
recovery operations may be further adapted according to the 
damping probability 7. 



I. Introduction 

In the most common treatments, quantum error correction 
(QEC) is developed for a very generic error model. An 
arbitrary error on a single qubit is correctable if both Pauli X 
and Z are correctable on that qubit; the continuum of quantum 
errors are thus reduced into a simple, discrete set. Using this 
approach, we may design error correction procedures that 
apply to a wide variety of quantum noise processes - the 
channel need only be well approximated by independent qubit 
errors. 

The general application of standard QEC comes with a 
price in efficiency. Quantum error correcting codes require 
a large number of redundant qubits; for short block lengths, 
generic codes are limited to low rates. While robust to arbitrary 
qubit errors, both error correction performance and efficiency 
can be improved by adapting the encoding and recovery 
operations to the physical noise process. Such adaptation is 
reasonable since, for any particular device, the noise will have 
a structure governed by the physical coupling of the system 
and the environment. Intuitively, we should be able to engineer 
improved error correction by careful adaptation of both the 
encoding and recovery operations. 

The concept of channel-adapted error correction is not new: 
early work labeled 'approximate' quantum error correction 
was presented in [1]. Much recent progress has been due 
to optimization efforts [2]-[7]. In each case, rather than 
correcting for arbitrary single qubit errors, the error recovery 
scheme was adapted to a model for the noise, with the goal 
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to maximize the fidelity of the operation. In [2], a semidefi- 
nite program (SDP) was used to maximize the entanglement 
fidelity, given a fixed encoding and channel model. In [5] 
and [6], encodings and decodings were iteratively improved 
using the performance criteria of ensemble average fidelity 
and entanglement fidelity, respectively. A sub-optimal method 
for minimum fidelity, using an SDP, was proposed in [7]. 
An analytical approach to channel-adapted recovery based on 
the pretty-good measurement and the average entanglement 
fidelity was derived in [8]. The main point of each scheme 
was to improve error corrective procedures by adapting to the 
physical noise process. 

The optimization efforts cited above detail mathematical 
and algorithmic tools with general application. That is to say, 
given any model for the noise process and an appropriately 
short code we can apply optimal [2], [5] and structured near- 
optimal [4] algorithms to provide channel-adapted encoding 
and recovery operations. 

It is important to note that the aforementioned tools are 
not, in themselves, complete solutions to the problem of 
channel-adapted QEC. When designing an error correction 
procedure, there is more to consider than whether an encoding 
or a recovery is physically legitimate. This motivated our 
exploration of near-optimal recovery operations [4], where 
we imposed a projective syndrome measurement constraint 
on recovery operations. Even given such a constraint, to 
implement channel-adapted QEC efficiently we need to design 
encoding and decoding procedures with sufficiently simple 
structure to allow efficient implementation. Furthermore, while 
the optimization routines focus on the entanglement fidelity 
and ensemble average fidelity due to their linearity, we should 
still like to understand the minimum fidelity, or worst case 
performance. 

To explore these issues in greater depth, we must consider 
channel-adapted QEC for a specific channel model. We exam- 
ine the amplitude damping channel, denoted £ a , given by the 
operation elements 
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Amplitude damping is a logical choice for several reasons. 
First of all, it has a useful physical interpretation: the param- 
eter 7 indicates the probability of decaying from state |1) to 
|0) (i.e. the probability of losing a photon). Second, amplitude 
damping cannot be written with scaled Pauli matrices as the 
operator elements; thus Theorem 1 from [3] does not apply 
and the optimal recovery operation does not have a near-trivial 
form. Finally, due to its structure, the amplitude damping 
channel can still be described with the stabilizer formalism, 
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TABLE I 

Optimal QER operator elements for the [4,1] code. Operators 

Rl AND R 2 CORRESPOND TO THE "NO DAMPINGS" TERM E^ 5 WHERE a 
AND DEPEND ON 7. R 3 - Rfj CORRECT FIRST ORDER DAMPINGS. 
R7 - RlO PARTIALLY CORRECT SOME SECOND ORDER DAMPINGS, 

THOUGH AS ONLY \0 L ) IS RETURNED IN THESE CASES SUPERPOSITION IS 
NOT PRESERVED. 



greatly aiding analysis. 

II. Qualitative analysis of channel-adapted QER 

FOR APPROXIMATE [4,1] CODE 

We begin with a qualitative understanding of the [4, 1] 
'approximate' code of [1] and its optimal channel-adapted 
recovery operation. The logical codewords are given by 
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Consider the optimal channel-adapted recovery for the [4,1] 
'approximate' code of [1] computed with the semidefinite 
program as presented in [2]. This is an example of a channel- 
adapted code, designed specifically for the amplitude damping 
channel rather than arbitrary qubit errors. Its initial publication 
demonstrated the utility of channel-adaptation (though without 
using such a term) for duplicating the performance of standard 
quantum codes with both a shorter block length and while 
achieving a higher rate. 

In [1], the authors proposed a recovery (decoding) circuit 
and demonstrated its strong performance in minimum fidelity. 
It is interesting to note that the recovery operation (described 
in quantum circuit form in Fig. 2 of [1]) is not a projective 
syndrome measurement followed by a unitary rotation as 
is standard for generic codes; instead it is a 7-dependent 
operation which includes a generalized (POVM) measurement. 
In contrast, the optimal recovery (in terms of entanglement 
fidelity) obtained via convex optimization does conform to 
such a structure. The optimal recovery operation is given in 
Table [I] We will analyze each of the operator elements in turn. 
For clarity of presentation, we begin with first and second 
order damping errors and then we turn our attention to the 
recovery from the 'no damping' term. 

A. Recovery for first and second order damping errors 

Neither Eq nor E\ in <[T} is a scaled unitary matrix, but 
we may understand the channel by considering E\ the 'error' 
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event. Let us denote a first order damping error as E\ , which 
consists of the qubit operator E\ on the k th qubit and the 
identity elsewhere. Consider now the effect of e[^ on the 
codewords of the [4, 1] code: 
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We see that the code subspace is perturbed onto an orthogonal 
subspace spanned by {|0111) , |0100)}. R3 projects onto this 
syndrome subspace and recovers appropriately into the logical 
codewords. Recovery operators R4, R$, and Rq similarly 
correct damping errors on the second, third, and fourth qubits. 
Notice that the first order damping errors move the information 
into mutually orthogonal subspaces. It is therefore not hard 
to see that the set of errors {/® 4 , E^}\ =1 satisfy the error 
correcting conditions for the [4, 1] code. (That the [4, 1] code 
satisfies the error correcting conditions for damping errors was 
pointed out in [9].) 

Consider now the subspace spanned by 
{J1010), |0101), |0110), |1001)}. By examining the logical 
codewords in Q and ([3j, we see that this subspace can only 
be reached by multiple damping errors. Unfortunately, in 
such a case we lose the logical superpositions as only |0l) is 
perturbed into this subspace. Consider, for example the two 
damping error Ef^f. We see that 
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While we cannot fully recover from such an error, we rec- 
ognize that these higher order errors occur with probability 
7 2 . Furthermore, we see that operator elements R? — Rxo 
do recover the |fjj,) portion of the input information. This 
contributes a small amount to the overall entanglement fidelity, 
though would obviously not help the minimum fidelity case. 
Indeed, R7 — Rxo do not contribute to maintaining the fidelity 
of an input state. 

We should also note that only a subset of all second order 
dampings are partially correctable as above. We reach the 
syndrome subspaces from Rt—Rxq only when a qubit from the 
first pair and a qubit from the second pair is damped, allowing 
the I Ox,) state to be recovered. If both the first and second 
qubits (or both the third and fourth qubits) are damped, the 
resulting states are no longer orthogonal to the code subspace. 
In fact, these are the only errors that will cause a logical bit 
flip, recovering |0l) as and vice versa. 

B. Recovery from the distortion of the 'no damping ' case 

We turn now to the recovery operators Rx and R 2 . Together 
these project onto the syndrome subspace with basis vectors 
{|0000) , 1 1111) , |1100) , |0011)} which includes the entire 
code subspace. We just saw that I® 4 together with single qubit 
dampings are correctable, but £f 4 does not have an operator 
element proportional to J® 4 . Instead, the 'no dampings' term 
is given by E® A which depends on the damping parameter 
7. Indeed, consider the effect of the no damping term on the 
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A standard recovery operation projects onto the code sub- 
space. Consider the effect of such a recovery on an arbitrary 
input state a |0l) +b The resulting (un-normalized) state 
is 
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The extra term \ distorts the state from the original input. 
While this distortion is small as 7 — > 0, both the original 
recovery operation proposed in [1] and the optimal recovery 
seek to reduce this distortion by use of a 7-dependent opera- 
tion. We analyze the optimal recovery operation for this term 
and compare its efficacy with the simpler projection. 

We see that R\ projects onto a perturbed version 
of the codespace with basis vectors {(a 1 0000) + 
/3 11111)), |0011) + ^ |1100))} where a and (3 are 
chosen to maximize the entanglement fidelity. We can use 
any of the numerical techniques of [2], [4] to compute 
good values for a and (3, but we would like an intuitive 
understanding as well, a and (3 (where \(3\ = \/l — |a| 2 ) 
adjust the syndrome measurement P\ so that it is no longer 
(0i| + (1l|, the projector onto the code subspace. 
If we choose them so that (0 L |Pi |0 L ) = (1 L \P 1 \1 L ) then 
we will perfectly recover the original state when syndrome 
Pi is detected for the no damping case. If syndrome Pi is 
detected, the no damping state will be distorted, but for small 
7, the second syndrome is a relatively rare occurrence. It 
could even be used as a classical indicator for a greater level 
of distortion. 

We can see in Fig. [TJ that the benefit of the optimal 
recovery operation is small, especially as 7 — > 0, though 
not negligible. Furthermore, the standard projection onto the 
code space is a simple operation while the optimal recovery 
is both 7-dependent and relatively complex to implement. 
For this reason, it is likely preferable to implement the more 
straightforward code projection, which still reaps most of the 
benefits of channel-adaptation. 

III. Amplitude damping errors in the stabilizer 

FORMALISM 

We turn our attention to the stabilizer formalism [9], [10], 
to demonstrate its utility for interpreting the amplitude damp- 
ing channel. In particular, understanding amplitude damping 
in terms of stabilizers allows a generalization of the [4, 1] 
code for higher rates. The stabilizer formalism provides an 
extremely useful and compact description for quantum error 
correcting codes. Code descriptions, syndrome measurements, 
and recovery operations can be understood by considering the 
n — k generators of an [n, k] stabilizer code. In standard prac- 
tice, Pauli group errors are considered and if {Xi,Yi, Z{yi =1 
errors can be corrected, we know we can correct an arbitrary 
error on one of the qubits since the Pauli operators are a basis 
for single qubit operators. 
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Fig. 1. Optimal vs. code projection recovery operations for the [4,1] code. 
We compare the entanglement fidelity for the optimal recovery operation 
and the recovery that includes a projection onto the code subspace. For 
comparison, we also include the original recovery operation proposed in [1] 
and the baseline performance of a single qubit. While the optimal recovery 
outperforms the code projector recovery, the performance gain is likely small 
compared to the cost of implementing the optimal. 



Let's consider the [4, 1] code in terms of its stabilizer group 

G = (XX XX, ZZII, 

II ZZ). We can choose the logical Pauli operators X — 
XXII and Z = ZIZI to specify the codewords in ^ and 
(pi, We saw in Sec. [n] that E^ damping errors together with 
™ 4 are correctable errors. Since each of these errors is a linear 
combination of Pauli group members: 
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(Xi + iYl), 



(11) 



we might presume that {J, Xi, Y{\f =1 are a set of correctable 
operations and the desired recovery follows the standard sta- 
bilizer syndrome measurement structure. This is not the case. 
Consider that the operator XiX 2 (or equivalently XXII) is in 
the normalizer N(G) of the code stabilizer, and thus {Xi, X2} 
are not a correctable set of errors. 

How, then, can the [4, 1] code correct errors of the form Xi+ 
iYil Instead of projecting onto the stabilizer subspaces and 
correcting Xi and Yi separately, we take advantage of the fact 
that the errors happen in superposition and project accordingly. 
As we saw, Xi + iYi and Xj + iYj project into orthogonal 
subspaces when i ^ j and we can recover accordingly. In 
fact, the correct syndrome structures can also be described in 
terms of stabilizers; understanding these syndromes enables 
design and analysis of other amplitude damping codes. 

Let G — (<7i, . . . , g n ~k) be the generators for an [n,k] 
stabilizer code. We wish to define the generators for the 
subspace resulting from a damping error Xi + iYi on the 
i th qubit. First, we should note that we can always write the 
generators of G so that at most one generator commutes with 
Xi and anti-commutes with 3^ (corresponding to a generator 
with an X on the i th qubit), at most one generator that anti- 
commutes with both Xi and Yi (corresponding to a generator 
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2 subspace 
-ZZI I 
I IZZ 
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4*' 1 subspace 
ZZI I 
-I IZZ 
1 1 IZ 



TABLE II 

Stabilizers for each of the damped subspaces of the [4, 1] code. 



with an Z on the i th qubit), and all other generators commute 
with both operators. Let |V) S C(G) be an arbitrary state in 
the subspace stabilized by G. If g £ G such that [<7, XJ = 
[g, Yi] = 0, then 



(Xi + iY,) |v) = (X + iY)g IV) = g(Xi + %%) 



(12) 



From this we see that the i th damped subspace is stabilized 
by the commuting generators of G. Now consider an element 
of G that anti-commutes with Xi and Y^ Then 

(Xi + iYi) |V) - [Xi + %Yi)g |V) = -g{Xi + iYi) |V) , (13) 

so — g is a stabilizer of the i th damped subspace. Finally, 
consider a g which commutes with Xi but anti-commutes with 

Y: 



(Xi + iY) |V) = (Xi + iYi)g |V) = g(Xi - iYi) 



(14) 



We see that neither g nor — g is a stabilizer for the subspace. 
It is, however, not hard to see that Z{ is a generator: 

Zi(Xi + iY) |V) = (iYi - i 2 Xi) |V) = (X t + iY,) |V) • (15) 

In this manner, given any code stabilizer G, we can construct 
the stabilizer for each of the damped subspaces. 

Consider now the stabilizer description of each of the 
damped subspaces for the [4, 1] code. These are given in Table 
[II] Recall that two stabilizer subspaces are orthogonal if and 
only if there is an element g that stabilizes one subspace while 
— g stabilizes the other. It is easy to see that each of these sub- 
spaces is orthogonal to the code subspace, as either —ZZII or 
—IIZZ is included. It is equally easy to see that the first and 
second subspaces are orthogonal to the third and fourth. To see 
that the first and second subspaces are orthogonal, note that 
—IZII stabilizes the first subspace, while IZI I stabilizes the 
second. Equivalently, —IIZI stabilizes the fourth subspace, 
thus making it orthogonal to the third. 

We can now understand the optimal recovery operation in 
terms of the code stabilizers. Consider measuring ZZII and 
IIZZ. If the result is (+1,+1) then we conclude that no 
damping has occurred and perform the non-stabilizer opera- 
tions of Ri and i?2 to minimize distortion. If we measure 
(—1,4-1) we know that either the first or the second qubit 
was damped. We can distinguish by measuring ZIII, with 
+ 1 indicating a damping on the first qubit and —la damping 
on the second. If our first syndrome is (+1,-1), we can 
distinguish between dampings on the third and fourth by 
measuring IIZI. If our first syndrome yields (—1,-1) we 



conclude that multiple dampings occurred. We could simply 
return an error, or we can do the partial corrections of R^—Riq 
by further measuring both ZIII and IIZI. It is worth 
pointing out a feature of the stabilizer analysis highlighted by 
this multiple dampings case. Each of the damping subspaces 
from Table [TT] has three stabilizers and thus encodes a 2 
dimensional subspace. Consider applying to the third 

damped subspace, equivalent to damping errors on qubits 1 
and 3. Note that there is no generator with an X in the first 
qubit; the resulting subspace is stabilized by 



-ZZII, -IIZZ, IIZI, ZIII). 



(16) 



As this has four independent generators, the resulting subspace 
has dimension 1. We saw this in the previous section, where 
for multiple dampings the recovery operation does not preserve 
logical superpositions but collapses to the |0l) state. 

Stabilizer descriptions for amplitude damping-adapted 
codes are quite advantageous. Just as in the case of standard 
quantum codes, the compact description facilitates analysis 
and aids design. While the recovery operations for the am- 
plitude damping codes are not quite as neatly described as 
the standard stabilizer recovery, the stabilizer formalism facil- 
itates the description. Furthermore, by considering stabilizer 
descriptions of the [4, 1] code and its recovery operation, we 
may design other channel-adapted amplitude damping codes. 
We will rely on stabilizers throughout the remainder of the 
paper. 

IV. Generalization of the [4,1] code for higher 
rates 

The stabilizer analysis for the [4, 1] code provides a ready 
means to generalize for higher rate code. Consider the two 
codes given in Table [TIT] (A). Each of these is an obvious 
extension of the [4, 1] code, but with a higher rate. Indeed the 
general structure can be extended as far as desired generating 
an [2(M + l),M] code for all positive integers M. We can 
thus generate a code with rate arbitrarily close to 1 /2. 



While the codes presented in Table III (A) have an obvious 
pattern related to the [4, 1] code, we will find it more conve- 
nient to consider the stabilizer in standard form as given in 
Table 



III 



(B). The standard form, including the choice of Xi 
and Zi, provides a systematic means to write the encoding 
circuit. The change is achieved through a reordering of the 
qubits which, due to the symmetry of the channel, has no 
effect on the error correction properties. 

Let's consider the form of the M + 2 stabilizer group 
generators. Just as with the [4, 1] code, the first generator has 
an X on every qubit. The physical qubits are grouped into 
M + 1 pairs; for each pair (i,j) there is a generator Z^Zj. 

The structure of the stabilizers makes it easy to see that 
j-j®2(M+i) )£! (fc)-}2(M+i) satisfy the error corre cting condi- 
tions for the [2(Af +1), M] code. To see this, we will show that 
the damped subspaces are mutually orthogonal, and orthogonal 
to the code subspace. Consider a damping on the i th qubit, 
where i and j are a pair. The resulting state is stabilized by Zi, 
—Z{Zj, and the remaining Z-pair generators. We will call this 
the i th damped subspace. This subspace is clearly orthogonal 
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Fig. 2. Circuit to encode the arbitrary state of M qubits given in qubits 
fcl • • • k M into 2(M + 1) physical qubits. This is the [2(M + 1), M] code 
in standard form. 



to the code subspace due to the presence of the —ZiZj 
stabilizer. For the same reason, the i th damped subspace is 
clearly orthogonal to the k th damped subspace for k ^ j. 
Finally, the i th and j th damped subspaces are orthogonal as 
we see that Z{ stabilizes the i th and — Zi stabilizes the j th . 

By writing the [2(M + 1), M] codes in the standard form, it 
is easy to generate an encoding circuit. The circuit to encode 
the arbitrary state \ip) in the M qubits k\ ■ ■ ■ kM is given in 
Fig. [2] The encoding circuit requires 3M+1 CNOT operations 
and one Hadamard gate. 

Let's write out the logical codewords of the [6, 2] code given 
the choice of Za in Table Hill 
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- |111111>) 
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—^=( 001001) H 


- |110110)) 
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- 1001111)). 
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Each codeword is the equal superposition of two basis states. 
We can see by inspection that the damped subspaces are 
mutually orthogonal: will eliminate one of the two basis 
states from each codeword and the resulting basis states do 
not overlap. 

A. Syndrome measurement 

We begin the recovery by first measuring the Z-pair sta- 
bilizers. A —1 result on the (i,j)-pair stabilizer indicates 
a damping of either the i th or j th qubit. This holds true 
even if multiple Z-pair stabilizers measure —1. Such a result 
indicates multiple damped qubits. Once we have identified the 
qubit pair, we perform an additional stabilizer measurement to 
determine which of the qubits was damped. As an example, if 
the (i, j)-pair was damped, we measure Zi, with a +1 result 
indicating a damping on the i th qubit and a —1 indicating a 
damping on the j th qubit. We perform this measurement for 
all pairs which measure — 1. 

If multiple stabilizers yield a —1 measurement then we have 
multiple damped qubits. As before, this reduces by half the 



XXXXXX 
(A) Z Z I I I I 

I I ZZ I I 
I I I I ZZ 



[6, 2] standard form 
XXXXXX 
ZZ I I I I 
I I Z I I z 
(B) IIIZZI 

Xi = III XX I 
X_ 2 = I IX I IX 
Zj. = ZI I I Z I 
Z 2 = ZI I I I z 



[8, 3] code 

xxxxxxxx 

Z Z I I I I I I 
I I ZZ I I I I 
I I I I ZZ I I 
I I I I I I z z 



[8, 3] standard form 
XXXXXXXX 
ZZ I I I I I I 
I I Z I I I I Z 
I I I Z I I Z I 
I I I I ZZ I I 
Xj = II I I XX I I 
X_ 2 = II I X I I X I 
X_ 3 = I IX I I I IX 
Zi = ZI I I I Z I I 
Z_ 2 = ZI I I I I Z I 
Z 3 = ZI I I I I I z 



TABLE III 

Stabilizers for [6, 2] and [8, 3] amplitude damping codes. In (A), 

THESE ARE WRITTEN IN A WAY TO ILLUSTRATE THE CONNECTION TO THE 
[4, 1] CODE. IN (B), WE PRESENT THE CODE IN THE STANDARD FORM, 
WHICH WE ACHIEVE MERELY BY SWAPPING THE CODE QUBITS AND 
CHOOSING THE LOGICAL OPERATORS SYSTEMATICALLY. THE STANDARD 
FORM PROVIDES A CONVENIENT DESCRIPTION FOR GENERATING 
QUANTUM CIRCUITS FOR ENCODING. 



dimension of the subspace and we cannot preserve all logical 
superpositions. For an example, examine the stabilizers for the 
[6,2] code when both the first and fifth qubits are damped: 



-ZZIIII, IIZIIZ, -IIIZZI, ZIIIII, IIIIZI). (21) 



This subspace has 5 stabilizers and thus has rank 2. Further- 
more, combining the last two stabilizers, we can see that 
ZIIIZI = Z\ stabilizes the subspace, indicating that the 
remaining logical information is spanned by {|01l) , |00i)}. 
In general, for a [2(M + 1), M] code, up to M + 1 dampings 
can be partially corrected as long as the dampings occur on 
distinct qubit pairs. If m is the number of damped qubits, then 
the resulting subspace has dimension 2 M+1 ~ m . 

If all Z-pair measurements for the [2(M + 1),M] code 
return +1, we determine that we are in the 'no dampings' 
syndrome and may perform some further operation to reduce 
distortion as much as possible. As in the example of the [4, 1] 
code in Sec. II-B we can choose to optimize this recovery 
a 7-dependent recovery or we can apply a stabilizer 



with 

projective measurement. In the former case, we may calculate 
an optimized recovery with a SDP or any of the near-optimal 
methods of [3], [4]. If we choose a stabilizer measurement, 
we simply measure the all-X generator X ■ ■ ■ X where a +1 
result is a projection onto the code subspace. A —1 result can 
be corrected by applying a Zi operation (in fact a Z on any 
one of the qubits will suffice). This can be seen by noting 
that the —X ■ ■ ■ X stabilizer changes the logical codewords 
by replacing the + with a — . 
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Fig. 3. Syndrome measurement circuits for the [2(M + 1) , Af] code. Circuit 
(A) measures each of the Z-pair stabilizers. If all of the measurements in (A) 
are +1, we are in the 'no damping' syndrome and we perform the syndrome 
measurement in (B). If the (i, j')-pair stabilizer measures —1, we perform the 
syndrome measurement in (C). 



B. Stabilizer syndrome recovery operations 

In the previous section, we described syndrome measure- 
ments to determine which qubits were damped. We also 
explained the extent to which multiple qubit dampings are 
correctable. We now present a straightforward set of Clifford 
group operations to recover from each syndrome. 

Consider a syndrome measurement in which we determine 
that m qubits ii, . . . ,i m were damped, where m < M + 1. 
We recover from this syndrome via the following three steps: 

1) Apply a Hadamard gate on the i\ qubit. 

2) With qubit ii as the control, apply a CNOT gate to every 
other qubit. 

3) Flip every damped qubit: X^ ■ ■ ■ Xi m . 

The procedure is illustrated as a quantum circuit for a two- 
damping syndrome and the [6, 2] code in Fig. [4] 

To see that this is the correct syndrome recovery for the 
[2(M + 1),M] code, we need to examine the effect of the 
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Fig. 4. Syndrome recovery circuit for the [6,2] code with the first and fifth 
qubits damped. 



three gate operations on the damped subspace stabilizers. In 
the syndrome where ii,...,i m are damped, we have three 
categories of generators for the resulting stabilizer group: — Z- 
pair stabilizers for the damped pairs, +Z-pair stabilizers for 
the non-damped pairs, and Z i± , . . . , Z im for each damped 
qubit. We need to see the effect of the recovery gate operations 
on each of these generators. Fortunately, we can demonstrate 
all of the relevant cases with the example of the [6, 2] code 
with the first and fifth qubits damped: 



-ZZIIII 
I IZIIZ 
-IIIZZI 
Z 1 1 1 1 1 
I IIIZI 



-XZIIII 
I IZIIZ 
-I IIZZI 
X 1 1 1 1 1 
I IIIZI 

YYXXXX 
I IZI IZ 
III ZZ I 

xxxxxx 

ZI I IZI 



CNOTi 's 



-YYXXXX 
I IZIIZ 
-III ZZ I 

xxxxxx 

Z I I IZI 

ZZI I I I 

I IZI IZ 

I I I ZZ I . (22) 

XXXXXX 

ZI I IZI 



The final two sets of stabilizers are equivalent since ZZIIII 
is the product of XXXXXX and YYXXXX. The first four 
generators of the resulting group are the code stabilizer. The 
last generator is Z\ which, as we saw before, indicates that 
the recovered information is spanned by {|00l) , |01l)} while 
the other two dimensions of information have been lost. 

While we have shown that the syndrome recovery operation 
returns the information to the code subspace, it remains to 
demonstrate that the information is correctly decoded. We 
can demonstrate this by considering the syndrome recovery 
operation on each of the Zi of the code. By showing that each 
of these is correctly preserved, we conclude that the syndrome 
recovery operation is correct. 

We have chosen the Z$ so that each has exactly two qubit 
locations with a Z while the rest are /. There are, therefore, 
five cases of interest. In case 1, neither of the damped qubits 
corresponds to a location with a Z. In case 2, the first damped 
qubit (ii) corresponds to a location with a Z. In case 3, one 
of the Z locations corresponds to a damped qubit, but it is 
not i\. In case 4, both of the Z locations correspond to a 
damped qubit, but neither is i%. Finally, case 5 is when both 
Z locations correspond to damped qubits and one is ii. 

Without loss of generality, we can see the effect of each case 
by considering an example using ZIIIZI and appropriately 
selected damped qubits. Consider case 1 where we can let 
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ii = 2: 



ZIIIZI 



In case 2, ii = 1: 
- ZIIIZI -^ H h 

_^X H — X{„ 



ZIIIZI — v.CNOT^ s ZIIIZI 

ZIIIZI. (23) 



-xiii z i ^ CN0T n' s -yxxxyx 



o 



(24) 



Notice that this last is equivalent to ZIIIZI as XXXXXX 

is in the stabilizer. 

In case 3, i\ = 2 while «2 = 5: 



ZIIIZI 



,x ( , ■••x i 



-ZIIIZI -> CN0T n' s -ZIIIZI 

ZIIIZI. (25) 

In case 4, let ii = 2,, 12 = 1, and 13 = 5^] 

ZIIIZI ZIIIZI ^ CN0T 'i ' s ZIIIZI 

_^.x 4l -x im ZIIIZI. (26) 

In case 5, ii = 1 and 12 = 5: 
ZIIIZI -^ H *i 

.X;, — Xi. 



XIII Z I ^ CN0T n' s yxxXFX 



(27) 



We see that in all cases, the recovery procedure correctly 
preserves the geometry of the encoded information, even in 
the case of multiple qubit dampings. It is worth emphasizing, 
however, that when multiple qubits are damped at least half 
of the information dimensions are lost. 

C. Performance comparison 

It is useful to compare the performance of each of the 
[2(M + 1),M] codes in terms of the damping parameter 7. 
Consider a comparison between the [4, 1] code and the [6,2] 
code. To make a valid comparison, we need to establish a 
common baseline. We do this by considering the encoding 
of two qubits with the [4, 1] code. For the completely mixed 
state p — 1/2, this is the equivalent of squaring the single 
qubit entanglement fidelity: 



F e (p®p,Ko£®TloE) = F e (p,IZo£f 



(28) 



This comparison is given in Fig. [5] (A). To compare multiple 
codes, it is more straightforward to normalize each to a single 
qubit baseline. This can be done by computing F^ 1 ^^ for an 
[n,k] code. The normalized performance for the [4,1], [6,2], 
[8,3] and [10,4] codes is given in Fig. |5](B). 

It is very interesting to note how comparably these codes 
maintain the fidelity even as the code rate increases. This 
is particularly striking when noting that each code can still 
perfectly correct only a single damping error. Thus, the [4, l]® 4 
can correct 4 dampings (as long as they occur on separate 
blocks) while the [10,4] code can only perfectly correct 1. Yet 
we see that the normalized performance is quite comparable. 

' While this contradicts our statement that the lowest numbered qubit would 
be ii, the assignment of i\ = 2 when the first qubit is also damped has no 
impact on the argument. 



[6.2] 

[22] 




Fig. 5. Performance omparison of generalized amplitude damping codes. 
In (A) we compare the [6, 2] code with the [4, 1] repeated twice. In (B), 
we compare the [4,1], [6,2], [8,3] and [10,4] codes. The entanglement 
fidelity has been normalized as 1/k where k is the number of encoded qubits. 
Notice that despite the increasing rates, the normalized entanglement fidelity 
maintains high performance. 



We take a closer look at the performance of the [8, 3] 
code in Fig. [6] We see that, while most of the entanglement 
fidelity is supplied by correcting no damping and E^ ' terms, 
a not insignificant performance benefit arises by partially 
correcting second order damping errors. In the case of the 
[4, 1] recovery, we concluded that such contributions improved 
the entanglement fidelity, but not the minimum fidelity as 
I lx) was never preserved by such a recovery. This is not the 
case for the higher rates. Two damping errors eliminate half 
of the logical space, but different combinations of damping 
errors will divide the logical space differently. For example, an 
damping error on the fifth and sixth qubits means the resulting 
space is stabilized by Z1Z2 thus eliminating logical states 
|01xl) and |10x^) (where x indicates either or 1). On the 
other hand, a damping on the fifth and seventh qubits results 
in a space stabilized by Z1Z3 eliminating logical states |0x1l) 
and \1x0l)- Thus, correcting second order damping errors still 
contributes to minimum fidelity performance. 

Given their identical rates, it is reasonable to compare the 
[8, 3] amplitude damping code presented here with the generic 
[8, 3] stabilizer code due to Gottesman [9]. The stabilizers for 



this code are presented in Table IV This code can correct 
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Complete [8,3] Recovery 
[3. J] Baseline 




Fig. 6. Fidelity contributions for each order error of the [8, 3] amplitude 
damping code. We see that the no damping, first, and second order recovery 
syndromes contribute to the entanglement fidelity of the recovery operation. 



Gottesman [8, 3] code 



XXXXXXXX 
ZZZZZZZZ 
I X I XY ZY Z 
I XZY I XZY 
IYXZXZ IY 



TABLE IV 

Stabilizers for the [8,3] code due to Gottesman [9]. 



an arbitrary single qubit error, and thus can correct all first 
order amplitude damping errors, as well as the less probable 
Z errors. These are corrected with 25 stabilizer syndrome 
measurements (Pauli operators on each of the 8 qubits as well 
as the identity). This leaves an additional 7 degrees of freedom 
to correct for higher order errors. While typically these are 
not specified, since we know the channel of interest is the 
amplitude damping channel, we can do a small amount of 
channel-adaptation by selecting appropriate recovery opera- 
tions for these syndromes. Since X and Y errors are the most 
common, we choose operators with 2 X's or 2 Y's (or one of 
each). 

The comparison between the rate 3/8 codes is given Fig. [7] 
Here we see that the channel-adapted [8,3] code outperforms 
the generic Gottesman code, but the effect is minor. The 
attention to higher order syndromes is seen to improve the 
performance of the [8, 3] code modestly. It should be pointed 
out that both recovery operations can be accomplished with 
Clifford group operations, and neither is dependent on 7. 

V. Linear codes for the amplitude damping 

CHANNEL 

The channel-adapted codes of the previous section have 
similar corrective properties to the [4,1] code: {I,e[ 1 '} are 
correctable errors while {Xi,Yi} are not. It is actually quite 
simple to design channel-adapted codes that correct both Xi 
and Yi errors and thus can correct {/, E\ ^} as well. Consider 
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Fig. 7. Comparison of the amplitude damping [8, 3] code and the generic 
rate [8, 3] code due to Gottesman. We include both the Gottesman recovery 
where no attention is paid to second order recoveries, as well as a recovery 
where second order syndromes are chosen to adapt to the amplitude damping 
channel. 



[7, 3] linear code 



I I I ZZZZ 
I ZZ I I ZZ 
Z I Z I Z I Z 
XXXXXXX 



TABLE V 

Amplitude damping channel-adapted [7, 3] linear code. Looking 

at the first three generators, this is clearly based on the 
classical Hamming code. The fourth generator differentiates 
between x and y syndromes. 



the [7, 3] code presented in Table [v] The first three stabilizers 
can be readily identified as the classical [7, 4] Hamming code 
parity check matrix (replacing with I and 1 with Z). They 
are also three of the six stabilizers for the Steane code. Measur- 
ing these three stabilizers, an Xi will result in a unique three 
bit measurement syndrome (Mi, Mj, M3), (In fact, a nice 
property of the Hamming code is that the syndrome, replacing 
+ 1 with and —1 with 1, is just the binary representation of i, 
the qubit that sustained the error.) Unfortunately, a Yi error will 
yield the same syndrome as Xi. We add the XXXXXXX 
generator to distinguish the two, resulting in 14 orthogonal 
error syndromes for the {Xi,Yi}J =1 . 

As in previous examples, we have a choice of recovery 
operations for the 'no dampings' syndrome. We can minimize 
the 'no damping' distortion as was done in previous cases 
by computing the optimal or structured near-optimal recov- 
ery within this subspace. This will result in a 7-dependent 
recovery operation. Alternatively, we can simply measure 
XXXXXXX with a +1 projecting onto the code subpsace 
and a —1 requiring a correction of Z,;. We compare these 
recovery operations in Fig. [9] 

We see in Fig. 10 that the [7, 3] code slightly outperforms 
the [8, 3] code of Sec. [fv] The [7, 3] code perfectly corrects 
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Fig. 8. Syndrome measurement circuit for the [7, 3] amplitude damping 
code. 
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Fig. 10. Comparison of the [7, 3] and [8, 4] qubit amplitude damping codes. 
We see that the [7, 3] performance is slightly better, despite the higher rate. 




Fig. 9. Optimal vs. code projection recovery operations for the [7,3] code. We 
compare the entanglement fidelity for the optimal recovery operation and the 
recovery that includes a projection onto the code subspace. For comparison, 
we also include the baseline performance of three unencoded qubits. While 
the optimal recovery outperforms the code projector recovery, the performance 
benefit is likely small compared to the cost of implementing the optimal. 



Fig. 11, Comparison of the [7,1] Steane code and the [7,3] amplitude 
damping code, normalized by 1 / k . We see that the [7, 3] performance is very 
similar to the EigQER optimized recovery for the Steane code. 



first order dampings and does not correct any second order 
dampings while the [8, 3] code partially corrects for higher 
order dampings. The performance advantage of the [7, 3] code 
arises from the decreased length: the probability of a higher 
order damping error decreases as only seven physical qubits 
are needed. 

Given its structure, it is logical to compare the [7, 3] 
amplitude damping code to the [7, 1] Steane code, as both are 
derived from the classical Hamming code. It was shown in 
[3], [4] that the Steane code is not particularly well adaptable 
to amplitude damping errors; despite its extra redundancy, 
the channel-adapted [5, 1] code significantly outperforms the 
channel-adapted Steane code. This is particularly unfortunate 
as the Steane code can be implemented with such efficiency, 



with particular value for fault tolerant quantum computing. 
The [7, 3] code provides a useful compromise position. 



We see in Fig. 1 1 the performance comparison for [7, 3] 
code and the [7, 1] code (with and without channel-adapted 
recovery). It comes as no surprise that the [7, 3] code out- 
performs the [7, 1] with standard stabilizer recovery: each 
perfectly corrects the first order damping errors, but the [7, 3] 
code has done so while preserving three times as much infor- 
mation. It is interesting to see how close the [7, 3] performance 
is to the channel-adapted [7,1]. It was shown in [3], [4] 
that the channel-adapted [7, 1] at least partially corrects some 
second-order damping errors; the [7, 3] code does not. This is 
mitigated by the higher rate of the [7, 3] code as again, three 
times as much information is preserved. 

The [7, 3] code is not the only high rate linear code for 
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amplitude damping errors. Consider any classical linear code 
that can correct 1 error for which codewords have even parity. 
We can convert this code to a quantum amplitude damping 
code in the same way as the [7, 3] code. If H is the parity 
check matrix for an [n, k] classical linear code, then each row 
can be made a quantum code stabilizer replacing l's with Z 
and 0's with /. To distinguish and Yj errors we include 
X®" as a generator. Since the classical code has even parity, 
we know that this generator commutes with the others. This 
construction yields a [n, k — 1] quantum amplitude damping 
code that corrects for single amplitude damping errors. 

The [7, 3] code we have presented here follows the structure 
proposed in [9] for amplitude damping codes; namely the code 
is a combination of an X-error correcting code and a Z-error 
detecting code. It is not immediately clear how to generalize to 
t error correcting linear codes. Instead of a single generator to 
distinguish Xi and Y{ errors, we require an extra t generators 
as we must distinguish X, and Yj for each corrected damping. 



Qubit 1 damped Qubits 2 & 3 damped 

-ZZI I I I I I I -ZZI I I I I I I 

IZZ I I I I I I IZZI I I I I I 

IIIZZIIII IIIZZIIII 

mi zz iii mi zz iii 

IIIIIIZZI IIIIIIZZI 

IIIIIIIZZ IIIIIIIZZ 

1 1 IXXXXXX 1 1 IXXXXXX 

ZI I I I I I I I -Z 1 1 I I I I I I 

Qubits 1 & 7 damped 
-ZZI 1 1 1 1 1 1 

IZZI 1 1 III 

1 1 IZZI 1 1 1 

1 1 1 IZZI 1 1 
-1 1 1 1 1 IZZI 

1 1 1 1 1 1 IZZ 

ZI 1 1 1 1 1 1 1 

1 1 1 1 1 IZI I 



TABLE VI 

Stabilizers for several damped subspace syndromes for the 
Shor code. 



VI. Amplitude damping errors and the Shor code 

We now turn our attention to the [9, 1] Shor code and 
its performance with a recovery operation channel-adapted to 
amplitude damping errors. The structured near-optimal results 
in [3], [4] showed that the Shor code provides remarkably good 
protection from the amplitude damping channel. Furthermore, 
the structured recovery operation is essentially optimal. Thus, 
in this case, the optimal channel-adapted recovery operation 
can be described as a projective syndrome measurement fol- 
lowed by a unitary operation. Given this intuitive structure, we 
can analyze the amplitude damping channel-adapted recovery 
operation. 

We first note that first order errors {E[ k ^} are perfectly 
correctable. This comes as no surprise, since the Shor code 
can correct an arbitrary single qubit operation. What may 
be surprising is that second order errors are also 

perfectly correctable. This was pointed out in [9] and can be 
seen through the same kind of stabilizer analysis of damped 
subspaces as we employed for the [2(M + 1), M] codes. 



We see in Table VI a few representative syndrome subspaces 
for damping errors on the Shor code. From these subspaces, we 
surmise that the first step in making a syndrome measurement 
is to measure the first 6 code stabilizers (each of which has 
a pair of Z's). Depending on those outcomes, we can make a 
further stabilizer measurement. 

As an example, consider when the first stabilizer returns a 
— 1 and the rest return +1. In that case, we can conclude that 
either the first qubit was damped, or both the second and third 
qubits were damped. These can be distinguished by measuring 
Z\ with a +1 indicating qubit one and a —1 indicating both 
qubits two and three. 

It is interesting to note that in this case, we will have only 
measured 7 stabilizers, and thus need one further measurement 
to achieve a 2 dimensional subspace. A natural choice would 
be to measure IIIXXXXXX, but this is an opportunity 
for a 7-dependent measurement instead. As before, such an 
operation can improve performance at the cost of circuit 
complexity. In most of this chapter, we have leaned toward 
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Fig. 12. Channel-adapted stabilizer recovery vs. 7-dependent recovery for 
the Shor code and the amplitude damping channel 



the simpler operation, concluding that 7-dependent operations 
provide some performance benefit but not enough to justify 
the added complexity. We will see that in the case of the Shor 
code, the performance gain may be sufficiently large to warrant 
a 7-dependent operation. 

Before this consideration, let's turn to another syndrome for 
multiple qubit dampings. We already examined an example 
where the second and third qubits are both damped. The Shor 
code is divided into three blocks of three qubits each; this case 
extends exactly to all circumstances when both damped qubits 

is an 



fall on the same block. The third example in Table VI 
example of two qubits damped from different blocks; in this 
case the first and fifth stabilizers are both measured to be — 1. 
While the most likely cause of this syndrome is a two-qubit 
damping, we can further measure Z^ and Z7 to correct for a 
three or four-qubit damping occurrence. 

The preceding discussion of stabilizer subspaces provides 
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two alternative recovery operations. We may begin with a 
projective syndrome measurement of the first 6 code gener- 
ators. At that point, we may either make a set of stabilizer 
measurements to project onto the damped subspaces, or we 
may make a 7-dependent syndrome recovery to minimize this 
distortion. It turns out that the best 7-dependent syndrome 
recovery has equivalent performance to the structured near- 
optimal recovery operation and is therefore essentially optimal. 
The stabilizer recovery, while simple to implement with Clif- 
ford group operations, has significantly weaker performance. 
We compare the two recovery operations for various values of 
7 in Fig. [12] 

How should we understand the extensive performance gains 
for the 7-dependent recovery? Both are consequences of the 
Eq distortion imparted onto the quantum state and the degrees 
of freedom in the code. The 7-dependent operation arises when 
we have a remaining degree of freedom after determining the 
syndrome. For the [2(M + 1),M] codes and the [7,3] code, 
we only have such freedom in the 'no damping' syndrome; 
in all of the damping syndromes, the syndrome measurement 
requires a full set of stabilizer measurements. We saw that for 
the Shor code first order dampings require only 7 stabilizer 
measurements to determine the syndrome, leaving one extra 
degree of freedom. We also have an extra degree of freedom 
when two qubits from the same block are damped. These 
constitute all of the first and some of the second order 
syndromes, each of which can be optimized to minimize Eq 
distortion. 

VII. Conclusion 

We have developed several quantum error correcting codes 
channel-adapted for the amplitude damping channel. All of 
the encodings can be compactly described in the stabilizer 
formalism. While optimized 7-dependent recovery operations 
are possible, a much simpler recovery operation using only sta- 
bilizer measurements and Clifford group operations achieves 
nearly equivalent performance. The channel-adapted codes 
have much higher rates (with short block lengths) than generic 
quantum codes. 

The creation of straightforward codes for the amplitude 
damping channel is a major step toward the still-open question 
of channel-adapted fault tolerant quantum computing. Intu- 
itively, channel-adaptation should be able to improve fault 
tolerant thresholds and reduce the necessary overhead for fault 
tolerance. Before such intuition is confirmed for the amplitude 
damping channel, several obstacles must be overcome. As an 
example, we must construct a universal set of operations on 
one or more of the channel-adapted codes presented here. 
Such challenges are acknowledged, but deferred for future 
consideration. 
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